Astronomy <fe Astrophysics manuscript no. 8740 


© ESO 2008 


February 2, 2008 





Letter to the Editor 



Inhibition of thermohaline mixing by a magnetic field in Ap star 
descendants: Implications for the Galactic evolution of 3 He 

Corinne Charbonnel 1,2 and Jean-Paul Zahn 3 

1 Geneva Observatory, University of Geneva, chemin des Maillettes 51, 1290 Sauverny, Switzerland 
e-mail: Corinne . CharbonnelOobs . unige . ch 

2 Laboratoire d'Astrophysique de Toulouse et Tarbes, CNRS UMR 5572, Universite Paul Sabatier Toulouse 3, 14, av. 
E.Belin, 31400 Toulouse, France 

3 LUTH, Observatoire de Paris, CNRS, Universite Paris-Diderot; Place Jules Janssen, 92195 Meudon, France 
e-mail: Jean-Paul.Zahn@obspm.fr 

Accepted for publication in A&A Letter, October 30, 2007 - To appear in Vol. 476, Issue 3 

ABSTRACT 

Aims. To reconcile the measurements of 3 He/H in Galactic HII regions with high values of 3 He in a couple of planetary 
nebulae, we propose that thermohaline mixing is inhibited by a fossil magnetic field in red giant stars that are descen- 
dants of Ap stars. 

Methods. We examine the effect of a magnetic field on the salt-finger instability, using a local analysis. 
Results. We obtain a threshold for the magnetic field of 10 4 — 10 5 Gauss, above which it inhibits thermohaline mixing 
in red giant stars located at or above the bump. Fields of that order are expected in the descendants of the Ap stars, 
taking into account the contraction of their core. 

Conclusions. We conclude that in a large fraction of the descendants of Ap stars thermohaline mixing does not occur. As 
a consequence these objects must produce 3 He as predicted by the standard theory of stellar evolution and as observed 
in the planetary nebulae NGC 3242 and J320. The relative number of such stars with respect to non-magnetic objects 
that undergo thermohaline mixing is consistent with the statistical constraint coming from observations of the carbon 
isotopic ratio in red giant stars. It also satisfies the Galactic requirements for the evolution of the 3 He abundance. 
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1. Introduction 

Numerous observations provide compelling evidence for a 
non-canonical mixing process which modifies the surface 
abundances of Li, C and N of low-mass red giant^j] when 
they reach the bump in the luminosity function (see refcr- 
\ ences in Charbonnel & Zahn 2007, hereafter CZ07). The 
quest for a mechanism which could lead to that extra- 
mixing between the hydrogen-burning shell and the con- 
vective envelope of red giant branch (RGB) stars has ini- 
tially been unsuccessful. For example rotational mixing, i.e. 
rotation-induced meridional circulation and shear turbu- 
lence (Zahn 1992; Maeder & Zahn 1998), is insufficient to 
explain the observed abundance pattern, as was shown by 
Charbonnel & Palacios (2004) and Palacios et al. (2006). 
But recently Eggleton et al. (2006) suggested a possible 
cause of such mixing, namely the molecular weight inversion 
created by the 3 He( 3 He,2p) 4 He reaction in the upper part 
of the hydrogen-burning shell. Based on 3D hydrodynamic 
simulations of a low-mass star at the RGB tip (Dearborn 
et al. 2006), they found that such a ^-profile is unstable 
and produces very efficient mixing. They claimed that this 
mixing was due to the well-known Raylcigh- Taylor insta- 
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1 i.e., stars with initial masses below ~ 2 - 2.5 Mq which evolve 
along the Red Giant Branch (RGB) to high luminosities until 
helium is ignited in their core under degenerate conditions. 



bility, which occurs in incompressible fluids when there is 
a density inversion. In stellar interiors, which are stratified 
due to their compressibility, a similar dynamical instability 
occurs when the Ledoux criterion for convection is satis- 
fied (Ledoux 1947), but it acts to render the temperature 
gradient adiabatic rather than to suppress the density in- 
version. Presumably it is this instability that Eggleton and 
colleagues observed with their 3D code, as attested by the 
high velocities they quote, because they started their sim- 
ulation with a unrealistic ID model which did not include 
any mixing, in which the /^-inversion had already attained 
its maximum value. 

In reality, the first instability to occur in a star, as the 
inverse /i-gradient gradually builds up, is the thermoha- 
line instability, as was pointed out by CZ07. This double- 
diffusive instability is observed in salted water in the form 
of elongated fingers, when the temperature is stably strat- 
ified, but salt is not, with fresh water at the bottom and 
salted at the top, the overall stratification being dynami- 
cally stable (Stern 1960). 

It was Ulrich (1972) who first noticed that the 
3 He( 3 He,2p) 4 He reaction would cause a ^-inversion, and 
he was the first to derive a prescription for the turbulent 
diffusivity produced by the thermohaline instability in stel- 
lar radiation zones. This prescription is based on a linear 
analysis, and it is certainly very crude, but it has the merit 
to exist. When it is applied to the ^-inversion layer in RGB 
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stars, with the shape factor recommended by Ulrich and 
sustained by laboratory experiments (Krishnamurti 2003), 
it yields a surface composition that is compatible with the 
observations of the carbon isotopic ratio as well as of the 
abundances of lithium, carbon and nitrogen in RGB stars 
(Fig. 3 of CZ07). 

Simultaneously thermohaline mixing leads to the de- 
struction of most of the 3 He produced during the star's 
lifetime (Fig. 4 of CZ07). It seems thus that we have fi- 
nally identified the mechanism responsible for abundance 
anomalies in RGB stars, which simultaneously accounts for 
the measurements of 3 He in Galactic HII regions (Balser et 
al. 1994, 1999a; Bania et al. 1997, 20020- 

A problem remains however, namely that all RGB stars 
brighter than the bump should undergo thermohaline mix- 
ing since they develop the same /i-invcrsion. As underlined 
by Balser et al. (2007), this was obviously not the case 
at least in two stars, i.e., in the RGB progenitors of the 
Planetary Nebulae NGC 3242 and J320. Slightly more mas- 
sive than the Sun, these PNe are indeed caught in the act 
of returning to the interstellar matter fresh elements syn- 
thesized in their womb among which 3 He with the amount 
predicted by classical stellar models (Balser et al. 1997, 
1999b, 2006; Galli et al. 1997). On the other hand the car- 
bon isotopic ratio data indicate that a few RGB stars escape 
that extra-mixing - about 4% according to Charbonnel & 
do Nascimento (1998). This is consistent with the predic- 
tions of the chemical evolution models, which suggest that 
fewer than 10% of all planetary nebulae enrich the ISM 
with 3 He at the level of classical stellar models (Galli et 
al. 1997; Tosi 1998, 2000; Palla et al. 2000, 2002; Romano et 
al. 2003), thus supporting the self-consistency of the extra- 
mixing mechanism. 

It is tempting to confront these numbers with the frac- 
tion of the Ap stars (suspected to all host magnetic fields) 
relative to all A stars, which is 5 — 10% as estimated by 
Wolff (1968). Since A-type stars are the progenitors of the 
more massive RGB stars, we are led to conjecture that in 
those giants that show no sign of extra-mixing, the thermo- 
haline instability is inhibited by a deeply buried magnetic 
field. 

Here we check that possibility by examining the effect of 
a magnetic field on the salt-finger instability, using a simple 
local analysis. 



2. Thermohaline instability in the presence of a 
magnetic field 

All perturbations (displacement magnetic field 6, pres- 
sure P' , temperature T", molecular weight p!) are expanded 
in Fourier modes 

exp[i(k x x + k y y + k z z) + st], (1) 

2 According to classical theory of stellar evolution which pre- 
dicts a large production of 3 He by low-mass stars, Galactic HII 
regions should be highly enriched in 3 He since they formed out 
of matter that has undergone 12 billion years of chemical evo- 
lution. However, their 3 He content is similar to that of the Sun, 
and thus no evidence for any enrichment of 3 He during the last 
4.5 Gyr was found. In addition, the "best" determination in a 
Galactic HII region (namely S209) by Bania et al. (2002) has 
yielded a 3 He abundance almost identical to the WMAP value, 
indicating that the stellar contribution to 3 He is very small. 



where we use local Cartesian coordinates, with the z-axis 
pointing in the vertical direction. 

We begin by deriving the buoyancy force, which is weak- 
ened through radiative damping; we neglect all other forms 
of dissipation that operate more slowly: Ohmic, molecular 
diffusion, viscosity. We split the buoyancy frequency in two 
parts, the first due to the thermal stratification and the 
second to the composition gradient: 



N 2 = Nf + Nf, 
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with the usual notations, and 
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We shall consider here only the case where the stratifica- 
tion is dynamically stable (N 2 > 0), i.e. when the Ledoux 
criterion for stability is satisfied: 
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The linearized heat equation may be written as 
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where k 2 = k 2 + k 2 + k 2 = k 2 + k 2 and k is the thermal diffu- 
sivity. (We simplify the Laplacian by assuming that k does 
not vary much over one vertical wavelength.) Likewise, the 
advection equation for the molecular weight perturbation 
takes the form 
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Thus the buoyancy force is given by 
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Turning next to the Lorcntz force, we perturb the mag- 
netic field Bq = (B x ,0,B z ) by the displacement £, and 
draw the field perturbation b from the induction equation: 



b = V x (£ x B ) 



(7) 



It is then straightforward to calculate the perturbation of 
the Lorentz force per unit volume: 



±- [(V x b) x B ] 



(8) 



where we neglect the variation of B over a perturbation 
wavelength. From here on we drop the index designating 
the basic field. 

It remains to implement the buoyancy and Lorcntz 
forces in the equation of motion, which we simplify by ne- 
glecting the effect of rotation. 



s 2 S + 



N 2 



1 + Kk 2 /s 



P' 

-ik h L. 

P 



(9) 
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Introducing the Alfvcn velocity a = B / ^Airp, we get 
P' 



ik x — + [s 2 + a 2 z {k 2 x + k 2 z )] 
P 

-\-kyd z (k x a z k z a x ) £y k x a x a z £ 2 = , 
P' 

iky 1- ky(k x a z — k z a x )(a z £ lX — a x £ z ) 

+ [s 2 + (k x a x + k z a z f + k 2 y (a 2 x + a 2 z )] f„ = , 

P' 2 2 

ik z — o, x a. z (k x + k z ) £ x kya x (k x a z k z a x ) £ y 

N 2 



(10) 



1 + KS 2 /a 



Nf, 



6 = 0, 



which we complete with the continuity equation, written in 
its simplest form, assuming that the Boussincsq approxi- 
mation is valid, i.e. that \k z \ ^> \d\np/dz\: 



k x £, x ~t~ ky^y ~t~ k z ^ z 



0. 



(11) 



The determinant of that fourth-order system is the disper- 
sion relation from which we derive the eigenvalues s, i.e. 
the growth-rate of the instability. 

In the non-dissipative case, this dispersion relation re- 
duces to (cf. Kumar et al. 1999) 



kf, 



(k-a) 2 )(s 2 + -±N 2 + (k-a) 



0. 



(12) 



The first set of (imaginary) roots corresponds to Alfvcn 
waves, and the second to internal gravity waves modified 
by the magnetic field. 

When thermal dissipation is included, these waves are 
damped, and a fifth root appears, which yields the growth- 
rate of the thermohaline instability in the presence of a 
magnetic field: 



-Kk 2 



N 2 (k h /k) 2 + (k ■ a) 2 
{N 2 + Nl){k h /kY + {k-af 



(13) 



Note that the magnetic field always plays a stabilizing role; 
it overcomes the destabilizing effect of the inverse fi gradi- 
ent when 



(k ■ a) 2 > -N 2 (k h /k) 



(14) 



We may assume that the horizontal size (A = 2ir/kh) of 
the 3 Hc fingers is much smaller than their vertical extent: 
kh ^ k z . Then it is mainly the horizontal component P>h of 
the field that prevents the thermohaline instability, when 



Bi> p —K\ 

7T 



(15) 



This criterion is sufficiently accurate for our purpose, al- 
though it has been established assuming that \N 2 \ does not 
vary much over a vertical wavelength of the perturbation, 
a condition which is not fulfilled in the real star. 

We are aware that the magnetic field itself may undergo 
instabilities (see Spruit 1999; Mcnou ct al. 2004) that are 
not captured by our crude analysis, since we ignored both 
rotation and the spatial variation of JBcEl- We have thus 



3 Busso et al. (2007) suggested that such instabilities might be 
the cause of extra-mixing. In that case however, the presence of 
deeply buried magnetic fields would have to be the rule, rather 
than the exception. We feel thus that the scarcity of magnetism 
observed in the A-type stars speaks in favor of the alternate 
scenario that we present here. 



Fig. 1. Critical value of the magnetic field which inhibits 
thermohaline mixing in a 1.5 M©, Z© model located at the 
bump. 



implicitly assumed that magnetic instabilities, if any, ap- 
peared before the inverse /x-gradient built up, and that they 
led to a stable field configuration. 

Figure 1 shows the critical value of P>h predicted by 
Eq. (|15p in the region of fi- inversion of a 1.5 M© model 
with solar metallicity. This model is located at the bump 
in the luminosity function, i.e., where thermohaline mixing 
should be fully efficient in the absence of magnetic field 
and where the C, N and Li abundance anomalies show up 
on the RGB (Charbonnel 1994). Here we assume that the 
finger width is 1/10 of the height of the /x-inversion region. 
According to our criterion, a magnetic field of £?/, f=a 10 5 
Gauss suppresses entirely the thermohaline mixing, while 
a field of 10 4 G is sufficient to inhibit the instability in the 
upper 1/3 of the \i inversion domain. We also computed 
(not shown here) the critical value of P>h at the bump in 
a 2 M© model at solar metallicity, and also in a 0.9 M© 
model with [Fe/H]=-1.8. Due to the fact that the structure 
of the H-burning region is very similar in all the cases, we 
obtain comparable values for B^ %t . 

Magnetic fields of that order are expected in the deep 
interior of the red giant descendants of the Ap stars. On the 
main sequence, these Ap stars possess surface fields of a few 
10 2 to about 3 x 10 4 G, which are believed to be of fossil 
origin; if so, the field is deeply rooted and its Ohmic decay 
time largely exceeds the main sequence lifetime. This sce- 
nario seems to be confirmed by a recent investigation which 
concludes that Ap stars whose mass is less than 3 M Q show 
no compelling evidence for field decrease during their main 
sequence life (Landstreet et al. 2007)0. The hypothesis of 
the fossil nature of magnetism in the Ap stars has recently 
been strengthened by the discovery that their progenitors 
are certainly the Herbig Ac/Be stars (Wade ct al. 2007a, 
b). 

After the main sequence phase, as the central region 
of the star contracts during evolution on the giant branch, 
the magnetic field may be enhanced by up to two orders 
of magnitude, due to flux conservation. Thus it can easily 
exceed the critical strength that ensures the suppression of 



4 The rapid field decline observed above this mass may be 
explained by the increased turbulent diffusion occurring in the 
larger convective core. 
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the thermohaline instability. We thus consider as plausible 
that those red giants which are prevented from undergoing 
the non-canonical mixing are the descendants of main se- 
quence Ap stars, because the magnetic field of fossil origin 
suppresses in them the thermohaline instability. If this sce- 
nario is correct, one should observe magnetic fields in RGB 
stars that behave classically, and in PNe such as NGC 3242 
and J320. 

3. Discussion 

The observational uncertainties in effective temperature, 
luminosity and metal content of Ap stars hamper a pre- 
cise determination of their mass (see Bagnulo et al. 2006). 
According to a survey covering the Ap/Bp stars located 
within 100 pc from the Sun (Power et al. 2006), magnetic 
stars range at solar mctallicity between 3.6 and 1.5 M© 
(with one exception at 6 M©). The lower limit coincides 
with the appearance of developed surface convection zones, 
which depends on the initial stellar metallicity. If this cor- 
relation held in the past independently of the bulk stel- 
lar metallicity, metal-poor magnetic stars could have ex- 
tended to lower initial mass: standard models predict in- 
deed that the convection zone of a 0.9 M© main sequence 
star with [Fe/H]=-2 has the same thickness as that of a 
1.5 M© main sequence star of solar metallicity (Talon & 
Charbonnel 2004). As a consequence, one can conjecture 
that below solar mctallicity, Ap star counterparts formed 
with initial masses much below 1.5 M©. This expected de- 
pendancy of the mass threshold for the magnetic inhibition 
of thermohaline mixing with metallicity is of direct rele- 
vance to the Galactic evolution of 3 He and will have to 
be taken into account in future chemical evolution models 
which consider metal-dependant yields. 

Galli et al. (1997) have determined a mass of 
1.2 ±0.2 M© for the progenitor of the planetary nebulae 
NGC 3242. Although such a determination is quite a feat 
and is undermined by several observational and theoretical 
difficulties, we find very encouraging that the initial mass 
derived for this object is so close to the lower mass limit for 
Ap stars at solar metallicity. We have to note indeed that 
the metallicity of NGC 3242 is half-solar (Barker 1985), so 
that the lower mass for Ap stars may already be shifted 
towards a lower value. 

It remains to be checked whether the magnetic stars are 
in sufficient number on the main sequence to account for 
the ~ 4% of red giants which arc spared from non-canonical 
mixing as estimated from the carbon isotopic ratio data 
(Charbonnel & Do Nascimento 1998). The classical value 
of the incidence of Ap/Bp stars, which is 5 — 10% based on 
a magnitude- limited sample (Wolff 1968), has been chal- 
lenged by the aforementioned survey performed by Power 
et al. (2006): in the solar neighborhood they find only 1.7% 
of magnetic stars. This value is close to that (i.e., ~ 3.5%) 
derived by North (1993) for Ap stars with masses between 
1.7 and 2.5 M© that are members of Galactic open clusters, 
and by Paunzen et al. (2006; ~ 2.2%) for chemically pecu- 
liar stars in young LMC open clusters. We note that the 
bulk incidence of Ap stars is very sensitive to the chosen 
value of the lower mass limit for these objects in a given 
sample (see Fig. 4 of Power et al. 2006 or Table VII of 
North 1993). In conclusion, we consider that the present 
statistics support our suggestion that thermohaline mixing 
is inhibited by a fossil magnetic field in those red giant stars 



that are the descendants of main sequence Ap stars. Within 
this framework we reconcile the long-standing problem of 
3 Hc overproduction on the Galactic timescale with the high 
values of 3 He observed in a couple of planetary nebulae. We 
encourage the observers to look for magnetic fields in "clas- 
sical" RGB and in NGC 3242 and J320 which we propose 
to call "thermohaline deviant stars" . 
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